function y = computeIsotropicHardeningModulus(element,alpha)
oY = element.TensileStrength;
K = element.IsotropicHardeningModulus;
y = oY+K*alpha;
end
